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Abstract. The nature of the probability distribution function of the local energy in 
the ±J Ising model has been investigated. At finite temperature, it has been derived 
that the probability distribution function must satisfy several relations at p = 1/2 
(p is the concentration of the ferromagnetic bond) and at Nishimori-line, respectively 
on any lattice in any dimension. They relate the probability distribution function 
corresponding to the local energy lower than — tanh(_ftT) with that corresponding to 
the local enegy greater than — tanh(-fT). [K is the inverse temperature.) The present 
results at Nishimori-line are, in a sense, generalization of Nishimori's result about 
the internal energy obtained by the local gauge transformation. Moreover, from the 
numerical calculation in the two-dimentional ±J Ising model, it is found that, in a 
certain temperature region, the probability distribution function of the local energy 
has several peaks which are related to the patterns of frustration around a bond of the 
lattice. 



PACS numbers: 75.50.Lk,64.60.Cn,05.50.+q 
1. Introduction 

To elucidate the nature of random spin systems, especially spin glass systems, has been 
a subject of a long-standing interest [1-21]. In the random spin systems, when we take 
one sample, namely, one bond configuration, the local energies of interacting bonds 
take various values, which change as the sample changes. There have been many works 
about the energy of spin glass systems, for example, about the ground state energy [7-9], 
energy barrier [10-12], energy landscape [13-15], low-energy excitation[16,17] for various 
spin glass systems. 

On the other hand, Nishimori[18] derived several rigorous results at a special line 
in the phase diagram of spin glass systems, which has now been called "Nishimori-line". 
The results were derived mainly by the use of the local gauge transformation, so that 
they hold for various spin glass models on any lattice in any dimension. The exact 
internal energy and the upper bound of the specific heat at Nishimori-line were derived. 
Moreover, a nature of the correlation function was derived, from which it follows that 
the ferromagnetic order parameter and the spin glass order parameter coincide with 
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each other at Nisimori-line, and the ferromagnetic phase boundary has some restriction 
in the phase diagram[18-20] . Furthermore, it was derived that the correlation-function 
distribution function must satisfy a certain relation at Nishimori-line[21]. 

In this paper, we investigate the nature of the probability distribution function of 
the local energy in the ± J Ising model. It has been derived that, at finite temperature, 
the probability distribution function mentioned above must satisfy several relations at 
p = 1/2 and at Nishimori-line, respectively, on any lattice in any dimension. They 
relate the probability distribution function corresponding to the local energy lower than 
— tanh(K) with that corresponding to the local enegy greater than — tanh(K). The 
present results at Nishimori-line are, in a sense, generalization of Nishimori's result about 
the internal energy [18], since the probability distribution function has more information 
than only the average value. 

Moreover, we have numerically calculated the probability distribution functions of 
the local energy at p — 1/2 and at Nishimori-line in the two-dimensional ± J Ising model, 
from which we have found that, in a certain temperature region, they have several peaks 
which are related to the patterns of frustration around a bond of the lattice. 



2. The probability distribution function of the local energy 

We consider the ± J Ising model, where the dimension of the lattice, the lattice structure 
and the range of the interactions may be arbitrary. The Hamiltonian is written as 
follows: 

H = -^TijSiSj, (1) 

(«) 

where Si = ±1, and the summation of (ij) runs over all the interacting pairs. Each 
is determined according to the following probability distribution: 

P(Tij) = p5(r t3 - 1) + (1 - p)5(r ij + 1). (2) 

In this paper, we put that J = 1 and &b — 1 (&b is the Boltzmann constant). 

Now, we denote the local energy of the interacting bond (ij) in a given bond 
configuration, {r}, as e^iK): 

eij(K) = - < TijSiSj > K , (3) 

where < ••• >k denotes the thermal average at temperature, T = 1/K. Then, the 
probability distribution function of the local energy, P e (x, K, K p ), can be written as 

P e (x,K,K p ) = [6(x-e ij (K))] Kp , (4) 

where [• • -]k p denotes the configurational average at the ferromagnetic bond 
concentration, p. (We define K p as exp(2K p ) = p/(l — p).) In this paper, we do 
not investigate the nature of P e (x ) K ) K p ) directly, but investigate that of a slightly 
different form, namely, 

P^(x, K, K p ) = [5(x - (cosh(2K) + smh(2K) ei ,(K)))] Kp , (5) 
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which is related to P e (x, K, K p ) as 

P e (x, K, K p ) = smh(2K)PP(cosh(2K) + sinh(2fT)x, K, K p ). (6) 

Here, it is noted that P e (x, K, K p ) may take non-zero value in the region, —1 < x < 1, so 
that Pj; 2 \x, K, K p ) may take non-zero value in the region, exp(— 2K) < x < exp(2K). 

3. The nature of the probability distribution function of the local energy 

In this section, we investigate the nature of the probabilty distribution function of the 
local energy, p( 2 \x, K, K p ), at p = 1/2 and at Nishimori-line, respectively. 
Firstly, we can derive the following identity: 

P^(x,K,K p ) 

= [(oosh(2tf p ) + sui H 2K p)eiJ (K p) )5(x - - {R) )] K , 

(7) 

(For the detailed derivation of equation (7), see the appendix.) 

Next, we integrate both terms of equation (7) from a to b (a and b are arbitrary 
real numbers, which satisfy exp(— 2K) < a < b < exp(2K).): namely, 

[ b p( 2 \x,K,K p )dx 

J a 

= j\cosh(2K p ) + sinH2K p )e tj (K p ))6(x - ^^^.^ M 

rl/a J. 1 dx 

= j i/b [(cosh(2Ay + *M2K, M K,M- - ^ (2 K ) + slnh(2JW ) l*F 

1/a 



f a [{cosh{2K p ) + sinh{2K p )e ij {K p ))5{x - (cosh(2fT) + smh{2K)e ij {K))] Kp dx, 

Jl/b 

(8) 



where we use the following property of 5-function to derive the last term: 

5(---)=x 2 5(x-a). (9) 

x a 

Equation (8) is the basic identity, which we use to derive the properties of the probability 
distribution function of the local energy. 

First, we investigate the case, p = 1/2. At p — 1/2, K p = 0, so that equation (8) 
becomes 

f b p( 2) (x,K,0)dx 

J a 

= [ 1/a [5{x - (cosh(2f0 + smh{2K)e ij {K)))} dx 

Jl/b 

p( 2 \x,K,0)dx (10) 
Changing the variable x into x^ 1 of rhs of equation (10), we also obtain 

[ b p( 2 \x,K,0)dx = [ b x- 2 p( 2 \x-\K,0)dx. (11) 

J a J a 
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For investigating the direct form of P( 2 \x, K, K p ), we define the following averaged 
distribution function, {P( 2 \x, K, K p } A : 

1 px-\-A 

{p( 2 \x, K, K P )} A = -J x Pj; 2 \x, K, K p )dx, (12) 

where A is an arbitrary finite value which satisfies x + A < exp(2K). Then, equation 
(11) can be written as: 

{P?\x,K,0)} A = {x- 2 PW(x-\K,0)} a . (13) 

We can take the value, A, as an arbitrarily small but finite value. In this sense, we 
can conclude that the value of P^(x, K, 0) coincides with that of x~ 2 P^ (x^ 1 , K, 0) at 
any finite temperature. Equations (10), (11) and (13) are the main results at p = 1/2. 
In a word, these equations are the conditions which relates the value of the probability 
distribution function at x with that at 1/x. As the temperature changes, the distribution 
function of the local energy itself may change, however, the distribution function must 
satisfy the above relations at any finite temperature on any lattice in any dimension. 
Using P e (x, K, K p ), equation (13) is written by the following form: 

{P e (x-(-tanh(K)),K,0)} A 

S 1 p( x-(-tanh(K)) 

~ l (l + smh(2K)(x- (-tanh(K)))) 2 e{ 1 + sinh(2A;)(x - (- tanh(X))) ' ' 



which becomes a rather complicated form, from which, however, we can see that equation 
(14) is a condition which relates the probability distribution function corresponding to 
the local energy lower than — tanh(fT) with that corresponding to the local energy 
greater than — tanh(i^). 

Next, we investigate the nature of the probability distribution function of the local 
energy at Nishimori-line. At Nishimori-line, K p = K, so that equation (8) becomes 

f b pW(x,K,K)dx 

J a 

= f 1/a [{cosh{2K) + smh{2K)e ij {K))5{x - (cosh(2fT) + smh{2K)e i:j {K))] K dx 

Jl/b 
rl/a 



ll/b 

Also, we obtain 

fb 



f a [x5{x - (cosh(2iT) + sinh{2K)e ij {K))} K dx 

Jl/b 

[ 1/a xP i2 \x,K,K)dx. (15) 

Jl/b 



f pW(x,K,K)dx= I x- 3 PP(x- 1 ,K,K)dx. (16) 

J a J a 

For the averaged distribution function, it can be derived that 

{PP(x,K,K)} A = {x-*p( 2 \x,K,K)} A . (17) 
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4. The property of the local energy at p — 1/2 and at Nishimori-line 

In this section, we derive several properties of the local energy at p = 1/2 and at 
Nishimori-line, using the results of the preceding section. 

For the configurational average of the local energy at Nishimori-line, from equation 
(15) , we obtain 

[cosh(2fT) + smh(2K)e ij (K)] K 

rexp(2K) 

= / xP( 2) (x,K,K)dx 

Jexp(-2K) 
r exp(2K) 

= / PP(x,K,K)dx = 1, (18) 

Jexp(-2K) 

from which, it is easily calculated that 

[ eij (K)] K = -tarih(K), (19) 

which was first derived by Nishimori[18]. At p — 1/2, by the similar procedure, it can 
be derived that 

[eij(K)] > -tanh(X). (20) 
Next, putting a = exp(— 2K) and b = 1 in equation (10) at p — 1/2, we obtain 

rl f exp{2K) 

/ pM(x,K,0)dx= / Pj: 2 \x,K,0)dx. (21) 

Jcxp(~2K) Jl 

Using P e (x, K, K p ), equation (21) can be written as 

r- tanh(JC) rl 

/ P e {x,K,0)dx = / P e {x,K,0)dx. (22) 

J-l J-tanh(K) 

On the other hand, at Nishimori-line, from equation (15), we obtain 

/ P e {x,K,K)dx> / P e {x,K,K)dx. (23) 

J-l J-tanh(X) 

The above results are intersting, since, at Nishimori-line, the configurational average 
of the local energy coincides with — tanh(X) at any temperature, while, at p — 1/2, 
the probability that the local energy, eij(K), takes the value smaller than — tanh(X) 
always coincides with one that eij(K) takes the value larger than — tanh(i^) at any 
temperature. Furthermore, these results hold on any lattice structure in any dimension. 



5. Numerical calculation of the probability distribution function 

Now, we show the examples how the above feature holds in the two-dimensional square- 
lattice ± J Ising model with only nearest neighbour interactions. By the transfer matrix 
method, we have calculated the values of eij(K) in L x L lattice (L = 11) for 10 7 
bond configurations, from which we have estimated {P( 2 \x, K, K p )}& at p = 1/2 and 
at Nishimori-line for A = (exp(2fT) — exp(— 2X))/1000. Figures 1 and 2 are the results 
of the above procedure at K = 0.75 at p = 1/2 and at Nisimori-line, respectively. We 
have confirmed that equations (13) and (17) definitely hold within the statistical errors. 
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In the figures, we can see that the averaged probability distribution function, 
{P( 2 \x, K, Kp)}^, has three peaks at p — 1/2, and one peak and two shoulders at 
Nishimori-line. In the two-dimensional square lattice ±J Ising model, the number of 
the patterns of furstration of two plaquettes around a certain bond of the lattice is 
three if we take into account only the number of frustration plaquettes, and each bond 
configuration belongs to one of the three patterns. We have confirmed that each peak or 
shoulder corresponds to the contribution from one of the three patterns of frustration. 
As the temperature decreases, each peak broadens, and at sufficient low temperature, it 
is found that {P( 2 \x, K, K P )} A has only one peak. More detailed numerical properties 
of the probability distribution function of the local energy is reported in a separated 
paper in the future. 

6. Conclusions 

We have investigated the nature of the probability distribution function of the local 
energy in Ising spin glass. 

We have derived that, in the ± J Ising model, the probability distribution function, 
P( 2 \x, K, K p )(P e (x, K, K p )), must satisfy several relations at finite temperature on any 
lattice in any dimension at p — 1/2 and at Nishimori-line, respectively. They relate the 
probability distribution function corresponding to the local energy lower than — tanh(i^) 
with that corresponding to the local enegy greater than — tanh(i^). The present results 
at Nishimori-line are generalization of Nishimori's result [18] about the internal energy, 
since the probability distribution function has more information than only the average 
value. 

Moreover, from the numerical calculation, we have found that, in a certain 
temperature region, the probability distribution function of the local enargy has several 
peaks which are related to the patterns of frustration around a certain bond in the 
lattice. More detailed numerical properties are reported in the near future. 
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Appendix A. Derivation of equation (7) 

In this appendix, we briefly explain the derivation of equation (7). 

First, we denote the partition function of the system in a given bond configuration, 
{r} as; 

Z(K,K') = 5> X P( E KT lm S l S m + K , T ij S i S j ), (A.l) 

{5} Im^ij 
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where we denote the temperature at bond (ij) separately. Also, we introduce the 
notation, [• • ] Kp K r. 

[' ' ~\k v ,k' v 

where A" B is the numbers of bonds, and we denote the ferromagnetic bond concentration 
of bond (ij) separately. Of course, if K p = K p , the notation mentioned above coincides 
with the standard one: 

[■■■}k p ,k p = [---]k p . (A.3) 
It is easily obtained that 

cosh(2X) + smh(2K) eij (K) = ( A - 4 ) 

Then, it yields that 

PP(x,K,K P ) = [5(x-(cosh(2K)+smh(2K) eij (K)))] Kp 

= [8(x - ^0y )U P ,K p (A.5) 

The last term of the above equation is invariant when we change the sign of K and K p at 
bond (ij) simultaneously, since it means that we just take the summation of Tij(— ±1) 
reversely; namely, 

PP(x,K,K p ) 

Z(K,K) 
= [8(x - z(Kj _ k) )U p ,-k p 

= [exp(-2K p Tij)S(x - z(k,—k) )\k p ,k p 

Z(K,K) 

= [(cosh(2Ag - ^(2K p ) nj )5(x - CQsh{2K) + ^ H2K) ^ {K) )U P 

= [(cosh(2Ag + s^(2K p )e i3 (K))5(x - - -L__ {R) )] Kp , 

(A.6) 

where we use the local gauge transformation to derive the last term. Thus, we obtain 
equation (7). 
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1.5 




Figure 1. The averaged probability distribution function, {Pe (x, K, 0)}a, at 
K = 0.75 &tp= 1/2. 
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Figure 2. The averaged probability distribution function, {Pe (x, K, K)}a, at 
K = 0.75 at Nishimori-linc. 



